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Abstract 

This paper presents new simple approaches for evaluating determinant and inverse of a matrix. The choice of pivot 
selection has been kept arbitrary thus they reduce the error while solving an ill conditioned system. Computation of 
determinant of a matrix has been made more efficient by saving unnecessary data storage and also by reducing the order of 
the matrix at each iteration, while dictionary notation [1] has been incorporated for computing the matrix inverse thereby 
saving unnecessary calculations. These algorithms are highly class room oriented, easy to use and implemented by students. 
By taking the advantage of flexibility in pivot selection, one may easily avoid development of the fractions by most. Unlike 
the matrix inversion method [2] and [3], the presented algorithms obviate the use of permutations and inverse 
permutations. 
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Introduction 

The term determinant is originally associated with the system of 
linear equations. It provides us with a foresight about the nature of 
solution of a given system of linear equations. 

Maclaurin [4] published the first result on determinant of 2 x 2 
and 3x3 system, which was generalized for n x n systems by 
Cramer [5]. Later, the well known Laplace expansion for 
evaluating determinant was proposed by Laplace, but he had 
then used the term 'resultant' instead of determinant. The first use 
of the term determinant in the modern context was done by 
Cauchy [6]. In 1866, Dodgson presented another method for 
finding the determinant of n x n systems which he named as the 
"Method of condensation" [7]. For larger systems the most 
preferred method for evaluating the determinant up till now is the 
Gaussian method [8] . It deals with the problem by converting the 
coefficient matrix into its equivalent upper/lower triangular form. 
The product of the pivot elements gives the determinant. 

The existence of matrix inverse also depends on its determinant. 
The method of finding the inverse by Gaussian method is 
discussed later in this paper. Ahmed & Khan [2] , and Khan, Shah, 
& Ahmad [3] proposed algorithms for the calculation of inverse of 
a matrix which are streamlined forms of the Gaussian method and 
also require permutation and inverse permutations. In this paper, 
we have presented two new algorithms related to the evaluation of 
determinant and inverse of a matrix. The first presented algorithm 
evaluates the determinant and is more efficient than Gaussian 
method as it reduces the order of the matrix at each iteration 
thereby saving unnecessary computations. In the second algo- 
rithm, we have presented another easy to mange way to calculate 
the inverse of the matrix by constructing the dictionary of the 



given system and thus excluding the need of permutations and 
inverse permutations. 

Determinant of A Matrix: A Brief Review 

Determinant of a square matrix A, denoted by det(A), is basically 
a real valued function. Because of its useful relationship with the 
matrix A and the solution to the system of equations of the form 
Ax = b, it becomes essential to have a knowledge about determi- 
nants while studying matrices. Evaluation of a determinant 
through its cofactor expansion (also known as Laplace expansion) 
is famous for lower order n x n matrices. 



Consider A = 



«11 «12 

a 2 \ a 22 
a ni a,a 



Let Mjj be the minor of entry a,y,(7e{l,2, ■ • • ,n}Je{l,2, ■ ■ ■ ,«}) 
which is the determinant of the sub-matrix obtained after deleting 
the i row and j th column of A. 

If / column of A is opted for cofactor expansion then, 



det (A) = ay Cy + a 2J C 2j + 
where, C,y is the cofactor of entry ay such that Cy = ( — Y) 1 My. 
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Similarly the cofactor expansion along the ('* row would be 
det (A) = a n C n + a i2 C a + ■■■ + a in C in 



For a matrix of order n, the evaluation of determinant by the 
above cofactor expansion requires computing n determinants of 
matrices of order («— 1). Therefore it can be implemented for 
finding determinant of a matrix of order 2 or 3 with ease but for 
higher orders it becomes a tedious job. To reduce the computa- 
tional effort usually the following three basic row operations of 
matrices are incorporated to evaluate the determinant [9], the 
method is known as evaluation of determinant by row reduction 
(also known as Gaussian method). 

Elementary row operations 

Let A be an nxn matrix, the following elementary row 
operations can be applied 

a) Multiply a row by a non-zero constant. 

b) Interchange two rows. 

c) Add a multiple of one row to another row. 



Effect of row operations on the value of determinant: [9] 

Let A be an n x n matrix then 

a) If B is the matrix that results when a single row or single 
column of A is multiplied by a scalar k, then def(B) = k det(A). 

b) If B is the matrix that results when a two rows or two columns 
of A are interchanged then det(B) = — det(A). 

c) If B is the matrix that results when a multiple of one row of A 
is added to another row or when a multiple of one column is 
added to another column, then det(B) = det(A). 



= — det 



■ — 2 det 



6 1 
-6 9 
1 5 

1 



property (b) 



3 2 



1 

3 -6 9 
0 1 5 



— 2 det 



0 -15 



= (-2)(— 15) det 



--(-2)(- 15) det 



property (a) 



1 

2 
15 

y 

0 1 5 



property (c) 



1 3 

0 1 
0 1 



1 3 
0 1 
0 1 



1 

2 

1 

~2 



1 

2 

1 

~2 
11 

y 



property (a) 



property (c) 



Lemma: [9] 

If A is an nxn triangular matrix (upper triangular, lower 
triangular or diagonal), then det(A) is the product of the entries on 
the main diagonal of the matrix; that is, det (A) = fliifl22 • • • a nn- 

Evaluation of determinant by row reduction 

The essence of the method is to transform the given matrix into 
its upper/lower triangular form by applying elementary row 
operations. The determinant can then be computed by incorpo- 
rating the properties defined above in 'Effect of row operations' 
and the lemma. 

Example. 



Consider the matrix A = 



1 5 

-6 9 
6 1 



det (A) = det 



0 1 5 

3 -6 9 
2 6 1 



= ( — 2)( — 15)(ll/2) lemma 3 



= 165 



The New Method 

Approach of row reduction may involve all the elementary row 
operations, as illustrated in above example. Here we are defining 
an operation (say pivot operation for evaluation of determinant) 
which consists of («— 1) applications of row operation (c) only, and 
can be employed to evaluate the determinants avoiding involve- 
ment of row operations (a) and (b). 

For example consider the following matrix, 



an ■■■ a y 



a\„ 



la„i 



PLOS ONE | www.plosone.org 



2 



February 2014 | Volume 9 | Issue 2 | e87219 



Determinant and Inverse of a Matrix 



Selecting a non-zero pivot element located at any arbitrary 
location say ay, and performing row operation (c) to make the 
remaining elements of the j th column zero we get 



an x ay 



a a 



a„i 



xa„ 



ay - 



c ay 



xa n 



«1„- 



:a y 



Step 5 : Reduce the order of A by removing p th row and k th 
column of A. Also set n: = n— 1 

Step 6: d : =(- \f +k x m x d. If «#0 go to step 2. 
Step 7: det(A) = d. Exit. 

Example 



Consider the matrix 



2 5 3 2 

4 10 1 7 

15 2 2 

2 12 1 



an 

a n xay 



ai„ — 



_ an 
a n \ x u n i 

a >j 



:ay 



a nn x a n 



Iteration 1: 

P. = {1,2,3,4}, here we take/>= 1 so, L: = {1,2,3,4}. Taking k= 1 



0 


-5 


5 


1 


2 


2 


-4 


-1 



Here we can see that the pivot row remains unchanged while 

the (p,k)' h element of the matrix, d pk = a p k ^ xa fc where a p j 

ay 

and a* are the elements corresponding to the pivot element in the 
pivot column and pivot row respectively. 

Cofactor expansion along the j' h column will give 

del (A) = OijCtj = ( - \) nj aijMjj 



det(A) = (-l) i+J aij det 



an 

a u x ay 

a v 



d\ n X Cl\j 



an 

a n \ x a n : 
a v 



The above procedure results in a computation of determinant of 
a matrix with reduced order of (n— 1). A repetition of the above 
procedure will eventually give a 2 x 2 determinant. For a »xn 
matrix there must be n non-zero pivot elements. The product of 
respective pivot elements gives the value of the determinant. 
However, if at any step there is no non-zero pivot exists, we can 
deduce that the determinant of the given matrix is 0. 

Problem 1 

Find determinant of the matrix A= [ay] of order n. 

Algorithm 1 

Step 1: Set d: = 1, 

Step 2: set R = { 1,2,3, ,«} 

Step 3: Select any peP such that L: = {k: a p k^=0, keP}. 
Step 4: If L = <p then d= 0 and go to step 7 
Otherwise m : =a p k peP, keL, and 



a J^Pl y ie p-{ p } and jeP-{k} 



d: =(-l) 1 + 1 x2 = 2 

Iteration 2: 

P. = {1,2,3}, here we take p=\ so, L: = {2,3}. Taking k=2 



5 13 
2 To 

-4 Z* 



d: =(-l) 1+2 x(-5)x2=10 

Iteration 3: 

P.= {1,2}, here we take p= 1 so, £ = {1,2}. Taking k= 1 



d : =(-l) 1 + 1 x(^)x 10 = 25 



Iteration 4: 

P. = {1}, here we take p= 1 so, Z: = { 1 }. Taking k= 1 

d: =(-l) 1 + 1 x(^)x25 = 12 
Hence determinant of given matrix is 12. 

Comparison with row reduction method 

The above example illustrates the efficiency of the algorithm in 
terms of memory storage and the number of elements computed. 
If a matrix of order 4 is solved by row reduction method the 
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number of elements computed will be 12, 6 and 2 in first, second 
and third iterations respectively. However we have computed 9, 4 
and 1 elements in the respective iterations. So the total number of 
element computations needed for row reduction method is 20, but 
our method needs only 14 element computations. Also, at each 
iteration the size of the matrix has been reduced. Row reduction 
method needs to store 1 6 elements for each iteration, hence on the 
whole method requires 48 elements to be stored in the memory, 
while in contrast our algorithm stores 16+9+4+1 = 30 elements, 
which is a noteworthy reduction in terms of storage requirement. 

A comparison of the number of elements computed and stored 
to evaluate the determinant of different orders by row reduction 
and our algorithm are shown in table 1. 

Inverse of A Matrix: A Brief Review 

Consider a matrix A= [ay] of order n. To evaluate the inverse of 
the matrix, say B, one must solve the following n system of 
equations, for X\ ,X2, ■ ■ ■ ,X n 



flllA'l a\2X2 
Cl2\X\ a 2 2X 2 



a n \X\ a„2X2 



Cl\ n X n 
A2//-V/ 



o 



system (1) 



From system (1) : x\ =b\\y\,X2 = 621^2, ■•■ ,X„ =b„\y„ 
From system (2) : x\ = buy\,x 2 = 622^2, • • • ,x„ =b„2y n 



From system (n) :x\ = b\ n y\ ,x 2 = b 2n y2, ■ ■ ■ ,X„ = b„„y„ 
bu b\2 ■■■ b\„ 



Then the matrix B = 



b%n 
b 1111 



is called the 



_b„\ b„ 2 
inverse of matrix A. 

This procedure could be compactly preformed by expressing 
the above system in the augmented form, 



Ax = I v 



(1) 

y„ ] T and I is 



where x=[xi x 2 ■■■ X n \'y=[y\ y% 
identity matrix of order n. 

If A is invertible, applying successive elementary row operations 
of Gaussian method yields, 



Dt = A~V 



a n \X\ a n 2X2 



Q.\ n X n 



0 

V2 



system (2) 



The New Approach Based on Dictionary Notation 

Now consider again the Equation (1) in expanded form, 



flu an 

fl21 A22 

a n i a„2 



fin 





~X\ 






X2 






_ "^n _ 





1 0 
0 1 

0 0 





'yi~ 




y% 




.y n _ 



Usually the following augmented matrix is used to solve above 
system, 



a\\X\ a\2X2 
a 2 \X\ a 22 x 2 



Cl\ n X n 
^2n 11 



d n \X\ a n 2X2 ' ' ' (lyi n X n 

The solutions of the above system are 



0 
0 

y n 



system (n) 



X\ X2 

flu an 

«21 «22 

a „i a„2 



n 


yi 


y 2 ■ 


■■ y n 


U\n 


1 


0 •• 


■ 0" 


O-ln 


0 


1 ■■ 


• 0 


Mini 


0 


0 •• 


■ 1 



Table 1. Comparison of row reduction method with algorithm 1 for evaluating determinant of matrix. 





Order 


Row reduction method 




Algorithm 1 






No. of Element Computations 


Storage Requirements (no. of 
elements) 


No. of Element Computations 


Storage Requirements (no. of 
elements) 


2 


2 


8 


1 


5 


4 


20 


48 


14 


30 


5 


40 


125 


30 


55 


7 


112 


343 


94 


140 


10 


330 


1000 


285 


385 


12 


572 


1728 


506 


650 
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Here, we are using the concept of dictionary notation developed 
by [1]. Now basic variables are the variables whose coefficients are 
in the form of any column of identity matrix, and a basis is 
collection of all basic variables. We can see that basis for the above 
matrix B= {y\,y2, • • • ,}>n} an d we may consider 
N = {x\,X2, ■ ■ ■ ,x n } as the non-basis of the current matrix. The 
objective is to convert the basic variables into non-basic variables 
and vice versa by using pivot operations. Using the dictionary 
concept defined by [1] we can remove the basic columns from the 
matrix and construct the following dictionary form with basis B 
and non-basis M (Note basic variables are shown in left most 
column and non basic variables in the top row of the dictionary). 



a p k 



1 

a P k 



Step 5: B: =B+{x k }-{y p } and N : = N-{x k } + {y p }. 
Update D(A) and go to step 2. 

Step 6: Inv(A) = [ay,Mi,jeH : XieB,yjeN\ Exit. 

Example 





X\ 


x 2 ■ 


Xn 




y\ 




an ■ 


■ a\„ 


Consider A 


yi 


an 


a 22 ■ 


■ a 2 „ 




y n 


_a n \ 


a„ 2 ■ 


Cl nn _ 





2 5 3 2 
5 10 1 7 

3 5 2 2 
2 12 1 



Pivot operation for evaluating inverse 

The following pivot operations [1] may be applied to enter x, 
into basis B and j>j into non-basis N, 

i) Divide the pivot row by pivot element and the pivot column 
by negative of the pivot element (except the pivot element). 

ii) The remaining («— l) 2 elements are determined by the 
formula as defined in the new method of determinant. 

iii) Reciprocate the pivot element. 

If the number of pivot elements is equal to the order of the 
matrix the resulting matrix gives the inverse otherwise we may 
conclude that the inverse does not exist. 

Problem 2 

Find the inverse of matrix A = [ay] of order n x n. 
Algorithm 2 

Step 1: Set H: = {1,2,3, ,«}, B: = {y h ) and N: = {x h ), 

MheH. Construct dictionary of the matrix A, i.e. -D(A). 
Step 2: Set P. = {p:y p eB} 
Step 3: If P= <p, go to step 6. 
Otherwise, L: = {t. a p k # 0, x^N) 
Step 4: If L= (p then inverse does not exist. Exit 
Otherwise, for any peP,keL 



m, : =- a ^\jieH-{p} 
a p k 



ay : = ay + a p j x m, V/e// — {p} ,VjeH — {k} 



a pj :=^VjeH-{k} 
a p k 





X\ 


x 2 


x 3 


x 4 


y\ 


'2 


5 


3 


2" 


yi 


5 


10 


1 


7 


y^ 


3 


5 


2 


2 


y4 


_2 


1 


2 


1 



D(A) = 



Here TV: = {x u X2,x 3 ,x 4 }, B: = {yuyi,y3,yi} 
Iteration 1: 

H: = {1,2,3,4}, R = {1,2,3,4}, taking p = 1 we get L: = {1,2,3,4}. 
Taking k = 1 



D(A): 



yi 
yi 
y* 



y\ x 2 x 3 

15 3 

2 2 2 

-2 0 -i 

1 5 1 

~2 2 2 

-1 -4 -] 



x 4 

1 
3 
1 

-1 



M = {yuX2,x } ,x 4 }, B-. = {xi ,y 2 ,yi,y4} 

Iteration 2: 

P. = {2,3,4}, taking p = 2 we get L: = {3,4}. Taking k= 3 



D{A): 





y\ 


x 2 


yi 


x 4 




- l 


5 


3 


19 " 




~Io 


2 


To 


To 


x\ 


2 


0 


l 


3 


x 3 


5 


~ 5 


~ 5 


73 


7 


5 


1 


13 


74 


~Io 

3 


2 

-4 


To 

l 


To 

8 




- ~5 




~ 5 


5- 



a,i : =m i 'iieH — {p} 



TV: = {yi,x 2 ,y 2 ,x 4 }, B: = {xi,X3,yi,y 4 } 
Iteration 3: 

P. = {3,4}, taking p= 3 we get L: = {2,4}. Taking k= 2 
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Table 2. Comparison of Gaussian method with the algorithm 2 for evaluating inverse of matrix. 



Order 


Gaussian Method 




Algorithm 2 






No. of Element Computations 


Storage Requirements (no. of 
elements) 


No. of Element Computations 


Storage Requirements (no. of 
elements) 


2 


12 


8 


8 


4 


4 


80 


32 


64 


16 


5 


150 


50 


125 


25 


7 


392 


98 


343 


49 


10 


1100 


200 


1000 


100 


12 


1872 


288 


1728 


144 
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D(A) = 





y\ 


>'3 


.1'2 


A'4 




- 3 


-1 

0 


1 


3 " 


X\ 


5 
2 


5 

1 


5 

3 


X3 


5 


~ 5 


~ 5 


x 2 


7 


2 


1 


13 




~25 
43 


5 
8 


25 

1 


25 
12 




-~25 


5 


~25 


25 - 



N = {yuyi,yz,M}, B-. = {xi,x3,x 2 ,y4} 

Iteration 4: 

P: = {4}, taking p = 4 we get L: = {4}. Taking k=4 



D(A) = 





y\ 


>'3 


J'2 


y4 




- n 




1 


5" 






-3 


4 


~4 


Xi 


7 




1 


5 


X3 


~4 


2 


~4 


4 


x 2 


19 


4 


1 


13 




12 


~ 3 


12 


12 


X4 












43 


10 


1 


25 




L 12 


3 


12 


12 - 



JVE = {yi,y3,yi,y4}, B-. = {xi,x 3 ,x 2 ,x 4 } 

Now place the elements with respect to indices of variables in i? 
and TV". 

For example Here H: = {1,2,3,4}, So xie.fi and y\eN, implies 
an =— . Also X\EB and 736^ implies an = — 3. Similarly placing 
the remaining elements we get 



D(A) = 





y\ 


yi 


73 


.V4 




- 11 


1 




5 " 




T 


4 


-3 


~4 


X\ 


19 


1 


4 


13 


x 2 


12 


12 


~ 3 


~T2 


X-i 


7 


1 


2 


5 




~4 


~4 




4 


X4 










43 


1 


10 


25 




-~ 12 


~ 12 


T 


12 - 



Comparison with Gaussian method 

The above example illustrates the efficiency of the algorithm in 
terms of memory storage and the number of elements computed. 
If a matrix of order 4 is solved by Gaussian method the number of 
elements computed will be 20 at each iteration. However our 
method requires computation of 1 6 elements at each iteration. So 
the total number of element computations needed for Gaussian 
method is 80, on the other hand for our method it is 64. Also the 
Gaussian method needs to store 32 elements for each iteration, 
hence on the whole Gaussian method requires 128 elements to be 
stored in the memory, while in contrast our algorithm stores 64 
elements, which is a noteworthy reduction in terms of storage 
requirement. 

A comparison of the number of elements computed and stored 
to evaluate the inverse by Gaussian method and our algorithm are 
shown in table 2. 

Applications 

Matrix determinant and inverse have applications in various 
fields like mathematics, economics, physics, biology etc. Solving 
different models like population growth involve the use of matrix 
determinant and inverse. Matrix inverse and determinant are also 
used in cryptography [10]. Linear transformations (rotation, 
reflection, translation etc.) involve the calculation of matrix 
inverse. Matrix inverse and determinant are also employed in 
operations research while solving linear programs, revised simplex 
method and markov chains. Determinants of order 3 are used to 
find area of triangles and for testing colinearity of points. Least 
square analysis of data requires the evaluation of matrix inverse 
[1 1]. /(-dimensional volume of parallelepiped in 5J m is determined 
by computing determinant [12]. 

Conclusion 

This paper presented easy algorithms for computations of 
determinant and inverse of a matrix. Since the order of the given 
matrix has been reduced at each step while calculating its 
determinant, the algorithm reduces the storage requirement (as 
exhibited in the example). The calculation of inverse has been 
done using the dictionary notation which obviates the use of 
permutations and makes it easier to cope with in class room 
teaching. Ill conditioned system can also be handled as the 
selection of pivots has been kept arbitrary, thus improving the 
numerical accuracy of the system. 
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